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FOREWORD 


The present report describes extensions to and generalizations of a radia- 
tion transport calculation procedure developed over a period of several years at 
first the Aerotherm Corporation, then the Aerotherm Division of the Acurex Cor- 
poration, Mountain View, California. The initial effort was done under Contract 
NAS 9-6719 for the NASA Manned Spacecraft Center, Structures and Mechanics Divi- 
sion. This effort included the development of the basic radiation properties 
and transport models. A subsequent effort was performed under Contract NAS1-9399 
for Langley Research Center, Applied Material and Physics Division. On this ef- 
fort, an equilibrium chemistry capability was incorporated into the procedure. 

The present effort was performed under Contract NASl-12,160 for the Langley Re- 
search Center, Space Systems Division. On this effort the radiation trans- 
port procedure was generalized to allow the equivalent width approximation, mo- 
lecular species properties associated with the Teflon and Silicon systems were 
added to the radiation model and the manner in which the far wings of the hydro- 
gen lines are calculated was modified. 
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SECTION 1 
INTRODUCTION 

I 

The ability to predict radiant energy transport is required in order to 
understand the heating phenomena associated with bodies entering planetary atmo- 
spheres at high velocities. The radiation transport is important both as an 
energy source in direct thermal contact and indirectly through the coupling be- 
tween it at the thermal boundary layer. Quantitative predictions of radiative 
transport under these conditions require frequency dependent properties of the 
radiating species, a frequency dependent transport model and a method for pre- 
dicting temperatures and species mole fractions, given the thermodynamic state 
conditions . 

In a prior study (Reference 1) , a procedure was developed for the pre- 
, diction of the appropriate thermodynamic equilibrium and radiation transport 
; events. This model has sufficient accuracy and requires modest computational 
efforts for applications where it is not coupled to a flow field procedure such 
as that of Sutton (Reference 2) . For the coupled application, the radiation 
transport has to be evaluated at space stations (diuny Lht; body) and at every 
i (or a good fraction thereof) iteration. Such requirements can cause the total 
; computational time to be measured in hours, even with modern computers. The 
| purpose of the present study is to modify the existing (Reference 1) radiation 
transport model to significantly reduce the computing efforts while maintaining 
acceptable accuracy. 

The radiation properties model is reviewed in Section 2; finite difference 
and equivalent width approaches for the line transport events are formulated in 
Section 3; a matrix of calculations is presented in Section 4; and concluding 
remarks are presented in Section 5. The calculations presented in Section 4 
illustrate the advantage in computing effort obtained and the trade-offs against 
accuracy which are required. 



SECTION 2 


RADIATION PROPERTIES 

The radiation properties model is based on that of Nicolet (Reference 1) . 
Two essential approximations are retained , the bandless molecular band system 
and the line group approximation for the atomic lines. The present model 
differs from that of Reference 1 in the treatment of the half widths of some 
of the hydrogen lines, as will be discussed subsequently. 

Molecular Band Model 

The bands within each band system are smeared according to the scheme 

% = / %dv/Av 

AV 

where the Av are selected such that the u v varies smoothly. The "bandless model" 
obtained in this fashion is felt to be satisfactory for radiation heating calcu- 
lations . 

Atomic and Ionic Line Model 

The line group approximation is employed to simplify the calculation. 

The frequency range of interest is divided into a number of frequency incre- 
ments - 15 to 20 or thereabout. These are not necessarily connected. Each 
one defines as a line group those lines which are centered within it. The line 
contribution at a frequency point within a frequency increment is obtained by 
considering only those lines within its group. This approximation is valid 
unless all three of the following conditions are satisfied: 

1. The path lengths are sufficiently long and the electron densities 
sufficiently high to lead to blanketing 

2. The line group boundaries are placed in between closely spaced lines 

3. The continuum intensities are not already black body 

In addition, the line group approximation becomes inaccurate for extremely long 
path lengths where the equivalent width of an isolated line can become greater 
than the width of the line group. This occurs frequently in hydrogen rich 
atmospheres where the wings of each of the four stronger hydrogen lines (Lyman 


3 



alpha and beta and Balmer alpha and beta) can significantly exceed the boundaries 
of its line group. This problem is solved by assigning the far wing contribu- 
tions to the continuum. The sketch illustrates the conditions where the line 
group approximation is valid (Figure la) , introduces small errors (Figure lb) 
and significant errors (Figure lc) . 



ACTUAL 

CONTINUUM 


Figure 1. Intensities in the Line Group Approximation 


Molecular Band Properties 

The y^ values are from the calculations of Woodward (Reference 3) , 
Biberman and Mnatsakanyan (Reference 4) and other reliable sources. 

Band systems are included for N 2 (l+,2+ and Birge-Hopf ield) , (Schumann- 
Runge) , NO (y, 3, 5 and e) , N^ + (1-), CN (Red and Violet ) , C 2 (Swan, 
Freymark and Mulliken) , H 2 (Lyman and Werner) and CO (4+) , all from the 
compilation of Nicolet (Reference 1). To this set, contributions have 
been added for the Swings band system, the SiO A^-X 1 ^ 4 * band system 
and the CF A 2 Z + -X 2 n band systems. 

The radiation properties model employs a Swings band cross 
section of 3x10 1 7 cm 2 which extends across the 3.0 to 3.5 ev region 
and which is temperature independent. This is consistent with the 
data of Brewer and Engelke (Reference 5) obtained in a furnace. Sufficient 
ficient information is not available to allow accurate extention of this 
data to high temperatures. The use of the low temperature cross sections 
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should be viewed as a method for obtaining estimates for the high temperature 
contribution. 

The contribution of the CF A 2 I + -X 2 n molecular band systems are 
from curve fits to radiation data generated with the Patch, Shackleford, 
and Penner "smeared line model" (Reference 6) , employing an f-number 
(of .0167) and Franck-Condon factors obtained from Wentink and Isaacson 
(Reference 7) . A sufficient set of molecular structure constants is 
available from the studies by Andrews and Barrow (Reference 8) , Porter, 

Mann, and Acquista (Reference 9) and Thrush and Zwolenik (Reference 10) . 

The contributions of the SiO A 1 n~X 1 E + molecular band system are 
also from curve fits to radiation properties generated with the smeared 
line model. The SiO f-number (of .13) is from the study of Smith and 
Liszt (Reference 11) , while the Franck-Condon factors are from the study 
of Liszt and Smith (Reference 12) . 

Atomic Line Properties 

With one exception, the hydrogenic lines considered individually 
employ line shapes from the studies of Griem (Reference 13) which consider 
electron and ion impact broadening. The exception is the half width 
relation for the high hydrogenic lines, derived originally as 


yJ, = const. N e log 10 


3 x 10 
n u N e 


f T 


+ n u 


.1/2 


which was found to be inapplicable to many of the conditions encountered 
during entry into the large planets. The problem was found in Griem' s use 
of a small argument expansion of an exponential integral during the 
original derivation. The elimination of the expansion yielded 


where 


y. = const. N 
x e 




V 


T 1 / 2 


<%<w 



( Y . ) 

v min' 



min 
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and 


y • 

^min 



which are suitable for numerical evaluation. 

Line properties are presented in Tables A.l, A. 2, and A. 3 in the ap- 
pendix for a radiation model suitable for entry calculations for most of the 
atmospheres suggested for Jupiter, Saturn and Uranus, and a second radia- 
tion model which is suitable for Venus is presented in Table A. 4. Some of 
the half-width data were obtained by estimate and are presented in paren- 
thesis . 
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SECTION 3 


TRANSPORT MODEL 


Formulation 

The basic equation governing the transfer of radiation through a medium 
in local thermodynamic equilibrium can be written as 


dl 

"ar = V B v 


V 


(i) 


where I is the spectral intensity, is the Planck function, S is the ray 
length, y^ is the absorption coefficient corrected for induced emission, viz. 

p v = M v ^ 1 " ex P( -hv / kT ) ] ( 2 ) 

and y^ is the ordinary absorption coefficient. 

In computing radiation fluxes across boundary and/or shock layers, it is 
convenient to make the tangent slab approximation. Thus, the properties along 
any ray can be related to those along the normal coordinate (y) by applying a 
cosine transformation, as shown in Figure 2. The resulting expressions for the 
optical depth, directional fluxes and total flux are well known and can be writ- 
ten in the form 


T 

v 



( 3 ) 


F >) 


-r 


V< )de t 


(4) 


F v (y) = 


/ 


V e v )de v 


+ 2J 


w 


& 


3 (t v ) 


(5) 
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OUTER BOUNDARY 



Figure 2. 


Radiating Layer Geometry 




( 6 ) 


/ » 

[F+ (y) - F~ (y) ] dv 

O 

where a diffuse wall has been assumed (reflectance r^ and emittance e™) and 
where and J™ are the black body emissive power and wall radiosity , respec- 
tively, and are defined as 


E 


v 


7TB 

v 


(7) 


J w = e w E (0) + r w F + (0) 

V V v v 1 V V 


( 8 ) 


Finally, the emissivities (e ) are 


+ 



1 - 



(9) 


where the t are dummy values of optical depth and the (r? (x) functions are 
exponential integrals of order n. 

The exponential approximation can be used to further simplify the equa- 
tions without an appreciable loss in accuracy for reentry problems. It is 
known (see Reference 14 for example) that when the approximation 



exp (-2x) 


( 10 ) 


is made, the transport solutions are exact in the optically thin limit, and of 
satisfactory accuracy in the optically thick limit. Thus, the emissitivi tes 
become 


V* 1 - exp(-2|? v - tj) (11) 

which are more convenient to work with than those given by Equation (5) . 
Emissivities written in this form have the additional advantage that by surpas- 
sing the factor of 2 in the exponential arguments and replacing E by B in 
Equations (4) and (5), the same forumlation can be used to calculate intensities. 

Additional simplification can be obtained for certain types of problems 
by employing a density stretched coordinate system. Introducing the 
transformation 
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£dn = p dy 


( 12 ) 


where 


C ¥■ S(y) 


yields 


% = C / K v dn < 13 > 

o 

for the optical depth. Here, the mass absorption coefficient 

K V E v V° (14) 

has been introduced. A final transformation of the form 

t = 2f 

introduces an optical depth variable which can be used directly as the exponen- 
tial argument in Equation (11) . 

Before evaluating the transport integrals, it is advantageous to separate 
them into line and continuum parts. This offers several advantages in that 
(1) optimum coordinates can be selected in frequency space for the finite dif- 
ference evaluation of the fluxes or (2) the transport integrals assume a form 
which is suitable for approximate evaluation by the equivalent width method 

and (3) an assessment of the importance of line effects is easily made from the 

c c • 

results of each given set of calculations. Substituting (or k for y^ xn 

~ c c 

Equations (2) to (5) yields corresponding values for and F . where the ± signs 
have been dropped from the latter for brevity. The line contributions to the 
flux are obtained by difference 


F 


L 

v 



(15) 


where the F values are evaluated using the total absorption coefficient (line 
plus continuum). Thus, the line contribution is treated as a correction to the 
continuum flux.* 

*It should be noted that the line contributions defined by Equation (15) can 
be negative for nonisothermal layers. 
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Finite Difference Evaluation of Line Fluxes 


In the evaluation of the transport integrals, it is assumed that the loga- 
rithm of the integrands can be represented as linear functions of the independent 
variables between adjacent spatial nodes. Thus, the mass absorption coefficients 
become (suppressing v subscripts for brevity) 


n - n. 


k ( TI ) = 



ni+i - ni 


(16) 


and the black body emissive powers can be written as similar functions of the 
optical depths. Substituting the interpolation functions into Equations (3) 
to (14) and performing the integrations yields the increment in the exponential 
argument 


AT i,i+l 


2 <iC(n i+ i 




K itl 

K , 

L 1 


£n 


li+i 

K . 

1 


(17) 


and recusion formulas for the directional fluxes 


< - •' 4,l ' 1+1 Ku * 


E .,"' 1 ' 1 * 1 - E . 


i+1 


1 + 


&n 


E . 
l 


lT i,i+l E i+1 


(18) 


, v iTi “ 1,i - E i-1 

F i * e ( r i-i + : r 


1 + 


£n 


At ... . E . , 

1-1,1 i-l 


(19) 


These equations are easily evaluated with two sweeps across the radiating layer 
provided that only one boundary is reflecting. The evaluation yields both 
directional fluxes at every node across the layer, consistent with the coupling 
requirements of most flow field procedures. Computationally, the use of the 
recursion formulas is much superior to a direct evaluation of Equations (4) 
and (5) and, under certain circumstances, comparable to the approximate equivalent 
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width approach, both of which require a sweep across the radiating layer for 
each pair of fluxes of interest. Moreover, the recursion formulas are suitably 
simple to be a point of departure for the development of influence coefficients 
(partial derivatives) required for some flow field coupling procedures. In 
view of these considerations, they were selected for numerical evaluation. 

Equations (17) to (19) are to be evaluated at a number of frequency points 
selected to adequately describe the variation of the spectrum. In the case of 
atomic lines, the number is typically 1500 frequency points (100 lines with 
15 points per line) , which is large enough to cause the evaluation of the line 
integrals to dominate the computational effort,. 

To obtain the maximum accuracy for a given computational effort, a fre- 
quency grid must be selected for each line and should be dependent upon the 
characteristics of the layer as well as the individual line. This is accomplished 
by introducing a parameter <J) which is characteristic of the width of the line 
(and will be discussed further in a subsequent paragraph) to stretch the coor- 
dinate system. The smallest increment and the distance to the most remote nodal 
point are defined by Equations (20) and (21), respectively, where v = |v - v | 

6v q = (20) 


v = ZJ 

max 2^ 


( 21 ) 


The quantities and z 2 are selected arbitrarily and are usually taken to be 
0.5 and 10 (respectively), on this order. The intermediate nodal points are 
established using a growth law, viz. 

<5 v . = ( 1 + z ) 6 v . n (22) 

1 D-l 

where the nodal spacing increases (increasing subscript j) with increasing 
distance from the center of the line. The growth factor z is determined implic- 
itly from the relation 


v 

max 




(23) 


where N is the number of increments to be used. Usually, the center of the line 
and about 5 to 7 additional points in each direction from it are sufficient. 
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The stretching parameter cf> is defined by Equation (24) 


4> - Y (?) yj 1 + t c (24) 

where the (half) half width y (y) is evaluated at a particular spatial location 
(usually as an edge condition) and t c is the optical depth of the entire layer 
at the frequency of the center of the line. It can be shown that (p has the fol- 
lowing properties ; 


lim <f) = y(y) 


T 

c 


0 


lim cj> * hv ~ 


(25) 


(26) 


y (y) / 


hv 



0 


That is, the <p is the distance in frequency space from the center of the line to 
the half intensity position when the line is weak, and it is approximately the 
distance to the frequency at which the layer has unit optical depth when the 
line is strong. At intermediate values of x , tj> should also be a reasonable 
approximation to the width of the line. Equations (20) to (26) specify the line 
frequency grid. 

Equivalent Width Evaluation of Line Fluxes 

Consider the absorption coefficient in a narrow spectral region 6v which 
contains an absorbing line. If 6v is sufficiently small, the continuum contribu- 
tion to the absorption coefficient y c can be approximated as 

u c f f(v) 

and the Planck function can be approximated as 

E v * f (v) 

with little loss in accuracy. If 6v is also sufficiently wide to contain the 
line, the incremental contribution to the flux* directed toward the wall is, i.e. 


*The remainder of this discussion will focus on the evaluation of the radiation 
flux, but is also applicable to the evaluation of the intensity. 
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AF + (0) 



6 v ° 


dt dv 


(27) 


can be written as 


AF + ( 0 ) = AF*(0) + AF+(0) 


where the continuum component has its usual definition 


AF+(0) = 6v / Ee -T c dt 
c o c 


(28) 


r y 

T = K 7 


M dy 

C g c y 


(29) 


and where y is the space coordinate measured from the point at which the flux is 
to be evaluated. The line component can be expressed as 


AF^ ( 0 ) 


-/ 


W(0) 


EdW 


(30) 


where the line equivalent width has been introduced 


W (0) = e 


/» - 

6 v 


) dv 


(31) 


where 


/• 


T i = K I ^ d y 


(32) 


but no additional approximations were imposed. The frequency integrations 
implicit in Equation (30) can be evaluated analytically for two interesting cases. 
Thus, Equation (30) reduces to 


AF 


y* 


i ' optically 
thin 


= / Ed[ exp(-t c ) j Sdy" 


(33) 
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for an isolated, optically thin line. Likewise, it reduces to 


Estrone, ' 2 /“I J/sW 1 (34) 

, . o o 

line 

for an isolated strong line.* 

In the Curtis-Godson approximation, the actual radiating layer is 
replaced by a homogenous layer which has approximately equivalent transport 
events. The homogenous layer is characterized by three parameters E^ , 
and $ h , where is the homogenous black body emissive power and the other 
two are defined as 




Vh 

2ity, 


= 


2try. 


(35) 


(36) 


where d is the mean line spacing and will be defined subsequently. It is 
reasonable to define the homogenous layer parameters by evaluating Equations 
(33) and (34) across a homogenous layer i.e. 


-t 




*«+l i-,*. 

Q | 


and 


which yield 


— T 


2E u e 

h 





strong line 


2 

IT 


AF* | optically thin 
AF+ | strong line 


2 


(37) 


(38) 


(39) 


*In this context, a strong line is defined as one in which the dominant 
radiative energy events occur in spectral regions far beyond the line 
half width, | v-v | >> Y . 
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( 40 ) 


Af| | optically thin 

$ h = — 

E h X h( e_T c/d) 

the desired results in terms of properties of the non-homogenous layer which 
can be evaluated directly. The flux emanating from a homogenous layer then 
becomes 

a K - Vh'w 1411 

where W h^ X h'^h^ the homo 9 eneous equivalent width. 

A variety of equations are available for predicting the equivalent width 
for a homogenous layer. Of these, possibly three are of interest to investi- 
gators studying radiation transport in the context of the entry environment. 
Thus, the equivalent width of an isolated line is given by the Ladenburg- 
Reiche formula 


W + = dB h L(X h ) 


(42) 


where LfX^) is the Ladenburg-Reiche function, which can be evaluated from well 
established exact or simple approximate analytic relations (References 15 and 16) 
or from tables. The equivalent width of an Elsasser band containing only strong 
lines is given by 


W 


+ 



(43) 


while the equivalent width of a band containing only weak lines is 

W + = d [1 - exp (“$ h x h )] 


(44) 


where a weak line is defined as being either optically thin or very strongly 
overlapped, after Plass (Reference 17) . The regions of validity of these three 
complementary relations are presented in Figure 3, which is consistent with the 
recommendations of Plass (Reference 18) . 

Valid estimates to the radiation flux increments from a real system can 
be made using Equations (41) to (44) provided that certain conditions are met. 
Thus, Equations (42) and (43) are valid only where the lines have a Lorentz 
shape. In addition. Equation (43) is valid only when the real line overlapping 


16 






events can be simulated through the us.e of the Elsasser model, a much more 
restrictive condition. However, it is felt that these conditions should be 
satisfied often enough for real systems to make the approach a useful one. 

The evaluation of the governing relations (Equations (33), (34), (39), 

(40), (41) and (42), (43), or (44)) is straight-forward except for the evalua- 

tion of AF^ | optically thin and AF* | strong line. As written,each quantity 
requires two integrations across the layer for each node in space where the 
flux is to be evaluated. This effort can be reduced for AF* | optically thin by 
transforming Equation (33) to a coordinate system measured from the wall, y, 
T C (Y) ' t c (y) , i.e. 


(AF «, I optically ) i 
thin 



y 

/ Sdy) 
o 


r 1 

J Sdy 



(45) 


which need be evaluated only once to yield the AF^ | optically thin at all the 
spatial nodes across the layer. The integrations in Equations (34) and (45) 
are readily evaluated by assuming that the logarithm of each integrand is a 
linear function of its independent variable between adjacent spatial nodal 
points . 

The governing relations for the radiative flux being transported away 
from the wall are similar to those presented here except for components due to 
wall reflection and emission which will be discussed subsequently. 

Equivalent Width of Wall Emission and Reflection 

The flux increment emanating from the wall includes components due to 
wall emission, continuum reflection and line reflection. Of these, the wall 
emission and continuum reflection components are treated in a straight-forward 
manner i.e. 


AF 


c 


w 


d(r w F^ 


w 


(46) 


for a wall characterized by a reflectivity r W , emissivity e W and black body 
emissive power E . The line reflection cannot be treated unambiguously under 
the existing assumptions, so it must be modeled. 
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Consider the flux increment incident upon the wall as being made up of 

a two step model (see Figure 4). The line component extends over a frequency 

increment W + where 
w 

W + £ d 
w 

and has a spectral flux of T - F+ . The continuum component consists of the 
usual spectral flux F + and extends W over the entire frequency range, d. With 
this boundary condition the optically thin flux increment becomes 


AF „ 


c y 

optically = / Ed 
thin J 


r y* " 

exp(-f c ) J Sdy " 


- F I- 
c 1 


+ r W (T - F + ) W + - min [w* 1,1 

w c W I W 11 

w 1 


(47) 


where 


1 1 = exp(-T c ) f^Sdy 

o 


Similarly, the strong line flux increment becomes 


AF , 


strong = 2 Ed 


/ i 

exp(-x )\ / Sydy " 

line o 

* o 


(48) 


-FI 
c 2 


+ r w (T - F + { W - mm 
v w c w 

w 




where 


= 2exp ( -T C ) 




SYdy ' 
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Figure 4. Flux Reflection Model 


As before, the optically thin and strong line flux increments can be employed 

to define the homogeneous layer flux. Four relations are possible for differing 

ratios of W + , I n and I 0 . 

w' 1 2 

Case Is I~ < W 4 *, I- < W + 

2 w 1 w 


X h 7T 


H 


10 


= 


h H 2 X h exp(-x c ) 


Case 2: I 2 < W+, ^ > W+ 


X. = - 

h 7T 


H 2 AF^ optically thin” 


(49) 


(50) 


(51) 


AF^ optically thin 
6 h = H 6 X h exp(-x c ) 


(52) 


Case 3: I 2 > wj, ^ 


x h = 7 


H H, 


H 2 AF^ strong line 


(b 3) 




H 2 X h ex P (_T c ) 


(54) 


Case 4: I _ > W + , I . > W + 

2 w 1 w 


x =1 
n 7T 


AF^ | optically thin 
AF„ [strong line 


(55) 


6 h = 


AF^ | optically thin 

H 2 X h eX P ( - T c ] 


(56) 
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where the coefficients are defined by 


H = AF 0 optically - H c 
thin 


(57) 


H_ = E - r W (T - F + ) - F 
2 h w c c 

w w 


(58) 


H_ = r W W + (T - F + ) exp (-T ) 
5 w w c ^ c 

w 


(59) 


H c = E, - F 
6 h c 


(60) 


H io = AF s> strong - h 5 

line 


(61) 


and are used in this form for all cases . 

Equations (47) to (61) can be evaluated once the spectral flux incident 
upon the wall and the frequency increment W* are evaluated. The usual equiva- 
lent width theory yields the product of these two quantities but not the quanti- 
ties themselves, so new information must be introduced. This is not a difficult 
task for a truly homogeneous layer as a simple relation of the form 


T - 
w 



, X 

(E - ) T“rV“ 

h w 1 X h 


(70) 


should be adequate. Note that Equation (70) is asymptotically correct for both 

large and small values of X, . At intermediate values of X, it should be of suf- 

n h 

ficient accuracy to meet the present requirements. Equation (70) was also se- 
lected for the general nonhomogeneous case. Its accuracy will depend on the 
effectiveness of the formulation in transforming a nonhomogeneous layer into a 
homogeneous one. Of course, the accuracy of the entire approach is strongly 
dependent upon this transformation, so that Equation (70) should not be viewed 
with any particular misgivings. 

Having defined the line component of the spectral flux by Equation (70) , 
the wall frequency increment is readily obtained from 


W 


+ 

w 


af; /(t, 

w 


w 


w 


( 71 ) 


which completes the two-step representation of the wall flux. 
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Selection of the Mean Line Spacing Parameter 

In an Elsassor band, the mean line spacing parameter d represents the 
distance (in ev) between the centers of adjacent lines. This parameter plus 
the mean line strength determines how much growth the line can experience 
before its far wings begin to overlap and interact with the far wings of the 
other lines in the band. Of course, an Elsassor band with its equally spaced 
lines of equal strength is only roughly similar to the real spectra which 
exists in the high temperature plasmas of interest here. Thus, the selection 
of d must be made to make the overlapping of the far wings of the real spectra 
look as much like as possible those of the lines in an Elsassor band. 

Consider the absorption coefficient at frequency due to the line 
centered at Then 


y 


i,k 


S . y. 

1 'i 


7T(V k - V.) 


(72) 


where the strong line approximation v, - v. >> y. has been invoked. At the 

K 1 1 

same frequency, the contributions from all the other lines in the band is 


y 


k 


£ 


S .y . 

1 JL 


7T (V. 


- V 


(73) 


If it is postulated that serious wing overlapping occurs when k' a 
significant frequency spacing parameter can be obtained trom Equations (72) 
and (73) as 


v i 


S i Y i 


z 

j^i 


-fill _ 

(V, - V.) 


(74) 


th 

which was adopted as one-half the mean line spacing for the i line, i.e. 


d. 


i 



Vi 


S .y . 
]] 



(75) 


This relation is based on strong line approximations and will not be valid 
in the weak line or optically thin limits; of course, the details of line 
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overlapping events are not important in these limits, so this should not 
be an important limitation on its applicability. 

The evaluation of with Equation (75) is straightforward for radiat- 
ing media which are uniform spatially. For nonuniform radiating media, either 
Equation (75) must be evaluated at conditions corresponding to those at a 
particular spatial point or some type of spatial averaging must be employed. 

In the present study, the first approach was adopted with the d^ being eval- 
uated at y. 

The Radiation Transport Program (RAD/EQUIL/1973) 

The entire solution procedure has been programmed for the UNIVAC 1108 
and the CDC 6600 and 7600 machines. A description of the programs has been 
presented elsewhere (Reference 19). Typically, the equivalent width approach 
allows a factor of four increase in computing speed. The accuracy of the 
equivalent width approach vis-a-vis the finite difference approach will be 
discussed in Section 4. 
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SECTION 4 


CALCULATIONS 


A matrix of theoretical predictions has been obtained and is presented 
in this section. The object was to obtain a thorough understanding of the 
strengths and weaknesses of the equivalent width approximation as compared to 
the finite difference approximation. To this end, calculations were made com- 
paring the fluxes of (1) an isolated line emanating from a homogeneous layer, 

(2) an isolated line plus a significant continuum emanating from a homogeneous 
layer, (3) an isolated hydrogen line plus a significant continuum emanating 
from a homogeneous layer, (4) an isolated line with a significant continuum and 
significant wall reflection, (5) a group of five identical lines in a line group 
emanating from a homogeneous layer, (6) a case including 20 line groups of 
varying features with the flux emanating from a nonhomogeneous layer, and (7) 
a general case with 20 line groups and a nonhomogeneous layer with varying 
elemental composition. 


Isolated Line in a Homogeneous Layer 

The predictions are presented in Table I for a radiating layer with 
optical properties of 


E = 
V 


Y = 
S = 


V 


C 


28,103 

0.0723 

8.42 x 10" 15 
0 


watts/ (cm 2 
ev 

ev - cm" 1 
cm“ 1 


ev) 


and where the width of the layer is treated as a parameter. For this problem, 
the equivalent-width approximation yields exact predictions. Thus, the - 5 per- 
cent comparison between the two sets of predictions should be viewed as a mea- 
sure of the truncation errors inherent in the finite difference approximation. 

In the predictions to be presented subsequently, a 5 percent comparison should 
be viewed as a meaningful limit, with better comparisons being fortuitous. 
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TABLE I 

HOMOGENEOUS LAYER 


( Layer 

Thickness 
| 6 
(cm) 


Line Flux Component 
at Wall (watts/cm 2 ) 
Finite 


Equivalent 

Width 


Difference 


Percent 

! 

Difference | 


j .001684 

.000057 

.000059 

3.4 

| .01684 

.00018 

.00018 

- 

I .1684 

.00060 

.00058 

3.3 

; 1.684 

.0019 

.0018 

5.3 

00 
I— l 

.0060 

.0059 

1.67 

! 168.4 

.019 

.019 

- 

! 336.9 

.027 

.026 

3.7 

842.2 

.041 

.040 

2.4 

Isolated Line Plus a Significant Continuum 


The predictions are presented in 

Table II for 

a radiating ! 

tical properties of 




E = 

V 

31,326 

watts/ (cm 2 - 

ev) 

Y = 4. 

79 x 10 -5 

ev 


S = 9, 

88 x lo" 3 

ev - cm " 1 


y c = 2. 

02 x 10" 3 

cm" 1 



and where the width of the layer is treated as a parameter. Observe that the 
line contributions are linear at the shortest path lengths, reach a maximum at 
intermediate path lengths, then decrease to negligible values at the longest 
path lengths. Thus, the continuum and line radiation do not interact signifi- 
cantly at the shortest path lengths, whereas the continuum totally dominates 
the radiative events of the layer at the longest path lengths. The ~ 6 percent 
comparison indicates that the presence of a significant continuum should cause 
a slight degradation in the accuracy of the equivalent-width approximation. 
Overall, the comparison is viewed as satisfactory. 
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TABLE II 

HOMOGENEOUS LAYER WITH A SIGNIFICANT CONTINUUM 


Layer 

Line Flux 

""*1 

Component 

Percent 

Thickness 

at Wall (watts/cm 2 ) 

Difference 

6 

Equivalent 

Finite 


(cm) 

Width 

Difference 



.00774 

.088 

- ■ ■ - i i 

.084 4.5 j 

.0774 

.82 

.78 

4.9 1 

.774 

4.8 

4.7 

i 

2.323 

8.1 

7.6 ' 

6.2 

7.74 

8.6 

8.2 

4.6 

23.23 

3.2 

3.1 

3.1 

! 77.4 

.026 

.025 

3.8 

774.0 

0.0 




Isolated Hydrogen Line 

The predictions are presented in Table III for a radiating layer with op- 
tical properties of 


E v = 28,103 

Y = 0.0723 

S = 5.16 

y = 2.31 x 10 -3 

c 


watts/(cm 2 - ev) 
ev 

ev - cm -1 


and where the width of the layer is treated as a parameter. Hydrogen lines, with 
their large wing spans and distinctive shapes, should be viewed as poor candi- 
dates for treatment with the equivalent width approximation. Indeed, the ~ 41 
percent comparison between the equivalent width and finite difference approxima- 
tion for the line contributions is poor. Consequently, the equivalent width 
approximation is not generally recommended for hydrogen lines. It should be 
observed, however, that the two approximations are in fair agreement ~ 20 percent 
as long as the lines are important contributors to the total radiation (line plus 
continuum) of the line group. This point is illustrated by the numbers in paren- 
theses in Table III which indicate fair accuracy ~ 14 percent for the calcula- 
tion of the total radiation assigned to the line group. Thus, one might argue 
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that decent estimates to even hydrogen line contributions can be obtained using 
the equivalent width approximation. 

TABLE III 


HOMOGENEOUS LAYER WITH A HYDROGEN 
LINE AND A SIGNIFICANT CONTINUUM 







Layer 

Line Flux 

Component 

Percent 

Thickness 

at Wall (watts/cm 2 ) 

Difference 

6 

Equivalent 

Finite 


(cm) 

Width 

Difference 


. 864 

130 

120 

7.6 

8.64 

1100 

960 

12.7 

00 

1700 

1500 

11.7 

129.5 

1100 

920 

16.3 (13.9) | 

172.7 

700 

540 

23. ( 6.3) | 

259.1 

230 

170 

26. ( .5)1 

518.2 

6 

4 

33. (.001); 

863.6 

.04 

.02 

41. — j 


Isolated Line with a Reflecting Wall 

The predictions are presented in Table IV for a radiating layer bounded 
on one side by a reflecting wall. The optical properties of the layer and wall 
are 


E = 

v 

Y = 
S = 


= 

r w = 


w 

£ 


2,310 

2.58 x 10~ 4 
17.86 

1.5 x 10~ 2 
0.6 
0.4 


watts/ (cm 2 
ev 

ev - cm" 1 
cm”" 1 


- ev) 


and where the width of the layer is treated as a parameter. Predictions are 
given at varous spatial stations (across the layer) for the line flux component 
directed toward the bow shock (away from the wall) . The ~ 10 percent comparison 
between the equivalent width and finite difference predictions represents some 
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additional decay in the accuracy of the equivalent width approximation. How- 
ever, the comparison is viewed as being satisfactory. 


Homogeneous Layer with Five Identical Lines 

The predictions are presented in Table V for a radiating layer with op- 
tical properties of 


Y 

S 


1,646 

1.48 x io -5 

5.4 x io -1 * 

4.5 x IO" 8 


watts/ (cm 2 - ev) 
ev 

ev - cm - 1 


where the five identical lines are equally spaced across a frequency increment 
of 0.1 ev and where the width of the layer is treated as a parameter. The 
line group responds as five individual strong lines at the shorter path lengths, 
as an overlapped line system at the larger path lengths. The - 8 percent maxi- 
mum difference occurs before overlapping effects become significant and can be 
assigned to inaccuracies which have been discussed previously. Thus, the ef- 
fects of line overlapping do not appear to cause significant decay in the ac- 
curacy provided that the Elsasser model is a reasonable representation of the 
actual line group. 


TABLE V 

HOMOGENEOUS LAYER WITH FIVE IDENTICAL LINES 


Layer 

Thickness 

6 

(cm) 

Line Flux Component 
at Wall (watts/cm 2 ) 

Percent 

Difference 

Equivalent 

Width 

Finite 

Difference 

! 

.1 

3.7 

3.6 

2.7 

.3 

6.4 

6.2 

3.1 

1.0 

12.0 

11.0 

8.3 

3.0 

20 . 0 

19.0 

5.0 

10.0 

33.0 

32.0 

3.0 
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THERMODYNAMIC PROPERTIES DISTRIBUTION 




Nonhomogeneous Layer with Real Line Groups 


The radiating layer is taken to have a width of 1.33 cm and thermodynamic 
properties as described in Table VI. There is assumed to exist a cold black wall 
bounding the layer on its low temperature side and a transparent shock wave bound- 
ary on the high temperature side. These are viewed as being typical of the in- 
viscid flow region about a blunt probe in the Venus entry environment. The op- 
tical properties are modeled in terms of 20 line groups which are described in 
the appendix. 

The important results are presented in Figures 5 and 6, where the fluxes 
toward the wall (Figure 5) and toward the shock wave (Figure 6) are presented as 
a function of spatial position for the total flux and a number of typical line 
groups. The comparisons are satisfactory, being within the 6 to 7 percent range 
for all the important line groups and for the totals. Moreover, the calculated 
flux at the wall is still more accurate, the comparison being within 2 to 3 per- 
cent. The equivalent width predictions tend to be consistently higher, so that 
one can assign some of the difference to truncation in the finite difference ap- 
proximation. Indeed, one might view these predictions as upper and lower bounds 
on the real transport events. 

The computing times for the UNIVAC 1108 are presented in Table VII. 


TABLE VII 
COMPUTING TIMES 


Computation 

Time 

(sec) 

Finite 

Difference 

Equivalent 
* Width 

Chemistry 

3.994 

3.994 

Continuum Radiation 

3.782 

3.782 

Line Radiation 

82.343 

7.082 

Total 

90.119 

14.858 


Thus, substantial savings in computational efforts can be obtained for a rela- 
tively small trade-off in accuracy. This should allow existing radiation-coupled 
procedures such as that of Sutton (Reference 1) to obtain comparable solutions 
for significantly reduced computational effort or, alternatively, allow consider- 
ation of additional physical events without significant additional computational 
penalty . 
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Line Group Contributions to Radiation Flux, Watts/cm 






Line Group Contributions to Radiation Flux, Watts/cnr 



Figure 6. Flux Directed From Wall 




Nonhomogeneous Layer with Varying Elemental Composition 

The radiating layer is taken to have a width of 1.452 cm and thermo- 
dynamic properties obtained from a thermodynamic equilibrium calculation and 
the temperature, pressures and elemental compositions given in Table VIII. 

There is assumed to exist a cold black wall on its low temperature side and a 
transparent shock wave boundary on the high temperature side. These are viewed 
as being typical of the stagnation region shock layer about an ablating probe 
in the Venus entry environment. The optical properties are modeled in terms 
of 20 line groups, including hydrogen lines, which are described in Reference 
19, under Sample Problem 5. 

The important results are presented in Figures 7 and 8, where the fluxes 
toward the wall (Figure 7) and toward the shock wave (Figure 8) are presented 
as a function of spatial position for the total flux and a number of typical 
line groups. Most of the line groups show excellent comparisons, especially 
at the lower frequencies (lower group numbers) . This suggests that the fin- 
ite difference and equivalent width could be combined to good effect for cer- 
tain types of problems. However, the line contributions to the flux differ 
by about 14 percent at the wall and about 15 percent in the interior, which 
is satisfactory for most flow field coupling problems. The continuum contri- 
butions to the flux values across the layer are comparable to the line contri- 
butions, so that the actual differences used in the flow field energy equation 
would be about one-half the values quoted above. 


35 




is to Radiation Flux, Watts/cm 



-°l 5fco 




,'95fcl 



SECTION 5 


CONCLUDING REMARKS 

A procedure has been developed for the evaluation of the line contribu- 
tions for an interesting class of radiation transport problems. The following 
points might be made in summary: 

• An equivalent width method has been formulated which considers 
continuum contributions and wall effects. 

• A number of calculations have been performed to build confidence 
in the method. 

• The comparisons with the finite difference procedure were all 
satisfactory with the possible exception of hydrogen lines. 

• A savings of considerable computational effort should be possible 
for many radiation-coupled flow field procedures. 
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APPENDIX A 


TABULATED PROPERTIES 


TABLE A. 1 


LINES AND LINE GROUPS FOR THE CHeH RADIATION MODEL 



Center 

Line Group j 

1 

~T 

Center 

Line Group 

Specie 

Frequency 

Boundaries 

i 

Specie j 

Frequency 

Boundaries 


(ev) 

(ev) 



(ev) 

(ev) 



<- .1 ^ 

, 

; 





H 

0.1662 



C 

1.487 



He 

0.1667 



He 

1.7024 


#6 

He 

0.2886 


Mi 


He 

1.7545 

| 


H 

0.3060 






-**1.7 9 j 

N 


H 

0.4722 

-*■ .55 



C 

1.814 





> 



He 

1.8562 


^ #7 

He 

0.5869 




H 

1.888 


He 

0.6022 




He 

2.1097 



He 

0.6493 






^2.4 J 


H 

0.660 


r #2 




~N. 


He 

0.6629 



He 

2.4714 



He 

0.6633 




He 

2.5185 



C 

0.6859 




H 

2.549 



C 

0.710 




He 

2.7722 



He 

0.7290 




He 

2.825 





o 

00 

+ 


H 

2.8559 


^ #8 






He 

2.9914 


He 

0.8218 




H 

3.0221 



C 

0.844 


y #3 


He 

3.0788 



C 

0.852 



He 

3.0918 





«- .94 

I 


He 

3.1266 





J 



He 

3.1877 



He 

0.9558 



He 

3 . 2453 



H 

0.9671 






^3.32^ 


He 

0.9691 







! 

He 

0.9695 


'>#4 


He 

3.3457 

1 

#9 

He 

0.9894 





-*-3.4 

C 

1.019 








He 

1.0356 




He 

3.4303 


! 

C 

1.0828 




He 

3.5955 


> #10 

H 

1.1333 






■*■3 . 65 „ 

i 

He 

1.1445 






i 


C 

1.163 

-**1.2 

; 








- 







C 

1.224 


i'#5 






C 

1.326 









TABLE A . 1 (Concluded) 


Specie 


Center 

Frequency 

(ev) 


Line Group 


Center 

Line Group 

Boundaries 

Specie 

Frequency 

Boundaries 


3.6952 

3.8886 

+ 4.0 

4.2090 

4.3816 

+ 4.475 

4.4851 

+ 4.8 

5.002 

+ 6.2 

6.424 

7.013 

7.078 

7.481 

7.717 

7.721 

7.947 

+ 8.0 ; 

8.03 

8.191 

8.203 

8.302 

8.368 

8.377 

8.433 

8.474 



+8 . 6 


9.137 

9.141 

9.332 

9.450 

9.611 

9.623 

9.697 

9.698 
9.709 
9.722 
9.797 
9.834 


10.196 

10.401 

10.405 

10.714 

10.873 

10.875 

10.887 

10.986 

11.061 


12.084 

12.181 


12.745 

13.052 

13.119 

13.2182 

13.601 


+ 11.2 


<- 11.2 


+ 13.8 J 








I 

iYMAN SERIES 


1-2 (L a ) 

1,215.67 

10.1968 

2 

0.4162 

1-3 (Lg) 

1,025.72 

12.0852 

2 

0.07910 

1-4 (L y ) 

972.54 

12.7460 

2 

0.02899 

1-5 (L 6 ) 

949.74 

13.0520 

2 

0.01394 

1-6 (L £ ) 

937.80 

13.2182 

2 

0.007799 

BALMER SERIES 

L ... 



2-3 (H a ) 

6,562.8 

1.8888 

8 

0.6407 

c a 

1 

CN 

4,861.32 

2.5499 

8 

0.1193 

2-5 (H y ) 

4,340.46 

2.8559 

8 

0.04467 

2-6 (H 6 ) 

4,101.73 

3.0221 

8 

0.02209 


3-4 (P a ) 
3-5 (Pg) 
3-6 (P Y ) 



18,751 

12,818 

10,938 



40,512. 

26,252. 

74,578. 


0.661 

18 

0.9671 

18 

1.1333 

18 

OTHERS 

0.306 

32 

0.4722 

32 

0.1662 

50 


0.8421 

0.1506 

0.05584 


1.038 

0.1793 

1.231 

















TABLE A. 3 


DATA FOR HELIUM LINES* 





Half 

Width 


Transition 

Line 

Center 

Frequency 

Per Free 
Electron at 
10,000 K 

Temperature 
Exponent ,* n, 
in Eq. 1 

Oscillator 

Strength 

f 



(ev) 

(ev-cm 3 ) 



Is3s-ls3p 

c : p) 

0.1667 



2 . 168 -2 ° 

-0.125 

0.1572 

Is3s-ls3p 

( 3 p) 

0.2886 

1.041 -20 

0.059 

0.672 

Is3p~ls4s 

( 3 s) 

0.5869 

- 2 0 

2.472 

0.292 

0.1087 

Is2s-ls2p 

(*p) 

0.6022 

5 . 924 -2 2 

0.386 

0.3764 

Is3p-ls4d 

( 3 D) 

0.6493 

l.iio" 19 

-0.274 

0.1617 

Is3d~ls4f 

(*F) 

0.6629 

d.o- 19 ) 

(-0.2) 

0.2525 

Is3d-ls4f 

( 3 F) 

0.6633 

(1. 0“ 1 9 ) 

(-0.2) 

0.1042 

Is3p-ls4d 

( 3 d) 

0.7290 

(1.110 -1 9 ) 

(-0.274) 

0.3615 

Is3s-ls4p 

( x p) 

0.8218 

9.541 -20 

-0.138 

0.035 

Is 3p-ls5d 

( ! D) 

0.9558 

(3.798~ 19 ) 

-0.325 

0.0347 

Is3d-ls5f 

('F) 

0.9691 

(3 . 0“ 1 9 ) 

(-0.2) 

0.0395 

Is3d-ls5f 

( 3 F) 

0.9695 

(3 . 0 - 1 9 ) 

(-0.2) 

0.0346 

Is3s-ls4p 

( 3 P) 

0.9894 

4 . 415~ 2 0 

0.031 

0.0321 

Is3p-ls5d 

( 3 D) 

1.0356 

(3 . 798“ 1 9 

(-0.2) 

0.0922 

Is2s-ls2p 

( 3 P) 

1.1445 

3 . 297~ 2 2 

0.560 

0.2695 

Is2p-ls3 s 

( ! S) 

1.7024 

(7.55“ ! ‘) 

0.331 

0.0480 

Is2p-ls3s 

( 3 S) 

1.7545 

5.539-“* ! 

■ 

0.331 | 

0.0347 

Is2p-ls3d 

( 3 d) 

1.8562 

- 2 0 

1.220 

-0.168 

0.711 

Is2p-ls3d 

( 3 D) 

2-1097 

6 S7«' 21 

n n 0 c 

V . 

A r~ r\ f\ 

kj • v\jy 

Is2s-ls3p 

( 3 P) 

2.4714 

„ „ 

- 2 0 

2.168 | 

-0.125 

0.1514 

Is2p-ls4d 

( X D) 

2.5185 

_ 1 9 

1.110 i 

-0.274 

0.122 

Is2p-ls4d 

( 3 D) 

2.7722 

(1.110- 19 ) i 

1 (-0.274) 

0.125 

Is2p-ls5d 

('d) 

2.825 

3 . 798~ 1 9 

-0.325 

0.0436 

Is2p-ls6d 

('d) 

2.9914 

(4.6“ 19 ) 

(-0.1) 

0.0213 

Is2p-ls5d 

( 3 d) 

3.0788 

(3 .798* 1 9 ) 

(-0.2) 

0.0474 

Is2p-ls7d 

( 3 d) 

3.0918 

(7.9 -20 ) 

(-0.2) 

0.0112 

Is2s-ls4p 

( 3 P) 

3.1266 

9.541- 20 

-0.138 

0.0507 

Is2s-ls3p 

( 3 P) 

3.1877 

1.041" 20 

0.059 

0.09478 

Is2p-ls6d 

( 3 D) 

3.2453 

(4.6 -19 ) 

(-0.2) 

0.0215 

Is2p-ls7d 

( 3 d) 

3.3457 

(7 . 9“ 1 9 ) 

(-0.2) 

0.0152 

Is2s-ls5p 

( ! P) 

3.4303 

2.506" 19 

-0.160 

0.0221 

Is2s-ls6p 

(*P) 

3.5955 

5.765“ 19 

! -0.173 

0.0128 

Is 2s-ls7p 

( 3 P) 

3.6952 

1.145“ 18 

| -0.173 

0.0066 

Is2s-ls4p 

( 3 P) 

3.8886 

4.415“ 20 

j 0.031 

0.0505 

Is2s~ls5p 

( 3 P) 

4.2090 

! 1 . 375“ 1 9 

! 0.005 

0.0293 

Is2s-ls6p 

( 3 P) 

4.3816 

3.253“ 1 9 

-0.021 

j 0.0169 

Is2s-ls7p 

( 3 P) 

4.4851 

1 

(7 . 9“ 1 9 ) 

(-0.2) 

0.0111 

1 


Numerical superscrips indicate multiplication by that power of 10 
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TABLE A. 4 


LINES AND LINE GROUPS FOR THE CON RADIATION MODEL 



Center 

Line Group 

Specie 

Frequency 

Boundaries 


(ev) 

(ev) 



«- .6 ^ 


0 

.685 



c 

.686 



N 

.689 


^#1 

C 

.710 



N 

.752 

— ~'\^- 
00 00 


c 

.844 



C 

.852 



N 

0 

.875 

.884 


>§2 

N 

.916 



N 

.930 

- .95 , 

- .96 S 


N 

.965 



0 

.991 



c 

1.019 



N 

C 

1.036 

1.083 


)>#3 

0 

1.098 


0 

1.132 



C 

1.163 

-1.2 , 


c 

1.224 


) 

N 

1.261 


1 

N 

1.319 


>#4 

C 

1.326 


0 

1.338 



N 

1.368 

-1.4 , 


N 

1.438 

> 

) 

0 

1.467 



c 

1.487 


>#5 

N 

1.553 



0 

j 

1.594 

-1.6 „ 




Center 

Line Group 

Specie 

Frequency 

Boundaries 


(ev) 

(ev) 




-1.62 ^ 


N 

1.663 



0 

1.767 



c 

1.814 


>#6 

N 

1.836 



0 

2.015 

-2.4 J 


N 

2.925 



0 

3.000 


>#7 

0 

3.167 


N 

3.472 

-3.5 J 
-3.4 "'j 


0 

3.711 

-4.0 J 

^#8 

C 

5.002 

■> 

-6.0 J 
-6.2 " 

>>#9 

C 

6.424 



C 

7.013 



c 

7.078 



N 

c 

7.111 

7.481 


>►#10 

c 

7.717 



c 

7.721 



c 

7.947 

^8. 0 ^ 


c 

8.030 

•N 


c 

8.191 



c 

8.203 



c 

8.302 



N 

c 

8.302 

8.368 


>#11 

c 

8.377 



c 

8.433 



c 

8.474 



N 

8.781 

4-9.0 y 













TABLE A. 4 (Continued) 



C 

N 

C 

N 

C 

N 

0 

C 

C 


C 

C 

C 

C 

N 

N 

0 

N 

C 

C 

N 


N 

N 

N 

N 

C 

N 

0 


Center 

Frequency 

(ev) 


9.139 

9.301 

9.332 

9.394 

9.450 

9.460 

9.501 

9.612 

9.697 


9.709 

9.722 

9.797 

9.834 

9.973 

10.102 

10.182 

10.332 

10.401 

10.405 

10.418 


10.493 
in hoc 

161619 

10.682 

10.714 

10.757 

10.761 


10.873 

10.875 

10.887 

10.927 




















TABLE A. 4 (Concluded) 



Center 

Line Group 

Specie 

Frequency 

Boundaries 

(ev) 

(ev) 



-<-13.8') 


N 

13.993 



N 

N 

14.160 

14.257 


>m 

N 

14.332 





-<-14.5 j 



NASA-Langley, 1975 


CR-2528 




